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f^i Abstract 



Since effective potentials derived from Feynman diagrams are naturally given in momentum 
space, we formulate the non-relativistic Coulomb problem entirely in momentum representation. 
We give momentum wave functions for all quantum numbers in one-dimensional integrals, even 
though they can be evaluated. Angular momentum decomposed Green's functions are then com- 
pactly represented. We apply this formalism to investigate the next to next leading order charm 
effects on IS bottomonium level shift. Our one insertion results are given completely in analytic 



form and numerically agree with previous results. Our two insertion results are also in agreement. 
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I. INTRODUCTION 

It is a common belief that new physics most likely will be due to new dynamics at a 
higher scale. Upon accepting this premise, we should expect that bigger effects will appear 
in processes in which known heavy particles, such as the bottom quark, the top quark and 
the Higgs boson actively participate. Such interplay between scales can be quite intricate, as 
exemplified in the p-parameter and in quite a few rare B decays. We would like on the one 
hand to use some of these processes to determine the masses of the known heavy particles 
as best as we can because the rates of many interesting processes depend on them strongly 
and on the other to lay the ground work to look for discrepancies as a byway to new physics. 
It is incumbent on us that for such processes the known dynamics must be well understood 
and that the so far uncalculable effects be well accounted for or at least be under control. 

One active area is the first few of bottomonium energy levels [1] and the production 
threshold of tt [2]. Here as the starting approximation, the heavy quarks move under a color 
Coulomb potential non-relativistically. Because retardation is small, one can add other 
effective interaction pieces as perturbations by integrating out the fast degrees of freedom. 
One can then extract out the bottom mass from the first couple of measured energy levels. 
For the top quark, it decays too fast to have real bound states, but these would-be bound 
states have strong effects on the shape of the production cross section near the thresholds. 

In the zeroth order approximation, such systems are just like a hydrogen atom. In most 
of the treatments of quarkonium systems, Coulomb-like wave functions in the coordinate 
space are used to evaluate matrix elements. We would like to ask whether an alternative 
may be also viable, for the reason that when one treats such systems in conjunction with 
quantum field theory, most if not all radiative correction calculations are done in momentum 
space via Feynman diagrams at some stage. In an effective Lagrangian approach, the Wilson 
coefficients are calculated in momentum space, but the operators are customarily converted 
back into coordinate space and their matrix elements are then obtained by averaging with 
spatial wave functions. It is therefore of some interest to see how to bypass the last step by 
following the momentum approach throughout. 

Needless to say, hydrogen wave functions in momentum representation have been re- 
peatedly used [3] for a long time. One difference we want to make here is that we have 



a convenient and compact representation for them in one-dimensional integrals, which we 
shall keep in such a form even though they can be evaluated into Gegenbauer polynomials. 
Matrix elements can then be evaluated mostly by the use of the method of residues and 
the result will collapse into only a small number of terms. We shall also give a very simple 
representation for the propagators. 

As we indicated earlier, we find it intriguing to explore the interplay of masses. Therefore, 
we shall use as an example to apply the momentum technique to investigate the charm mass 
effects [4 — 8] on the bottomonium level shifts to the next-next leading order. One notices 
that the binding momentum « |a s mb p» 2 GeV is comparable to m c ~ 1.25 GeV. One must 
take the charm loop as a whole as a term in the potential, which becomes highly non-trivial 
in the coordinate space. 

There are two sets of diagrams. The first set is due to the insertions of the lowest order 
charm loop once or twice, or the limiting zero mass particle effects and a charm loop each 
once on the same gluon line ( Fig.(l)). Also, there are diagrams due to the fourth order 
charm loop. (We neglect the charm loop effects on the vertex.) We have been able to 
evaluate all of them analytically into simple functions. They account for 98% of the charm 
effects for the IS level shift. The second set is due to double insertions of the above potential 
on two separate gluon lines ( Fig. (2) ) and so far can be given only in terms of some simple 
integrals. The results are in agreement with those given in p|. Furthermore, because we 
have the explicit functions for the first set at hand, we can follow the relevant branch of 
the functions and continue them to cover the toponium would-be bound state energy level 
shifts, here due to the bottom quark mass effects. We hope that this example is convincing 
enough to illustrate that the same technique can be used to cover exotic particle effects, 
when called upon. 



In an article |9J co-authored by one of us, we showed how the technique presented here 
could be extended to yield results at least for the spherically symmetrical states of the 
hydrogen atom and gave a bound on the non- commutative scale in a certain extension of 
the non-relativistic kinematics [10], by using the highly accurate 1S-2S energy difference. 



The plan of this paper is as follows: in t 



wave functions. An operator algebra 
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re next section, we shall solve for the momentum 
will be used to normalize them. These wave 







FIG. 1: We show schematically the next leading order and next-to-next leading order charm mass 
effects in single insertions which alter the bottomonium energy levels. The loops are composed of a 
charm and anti-charm pair, the shaded ellipses represent cumulative second order zero mass quark 
and gluon effects, and the double lines depict either a b or b quark. 
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FIG. 2: We show schematically the next-to-next leading order charm effects in double insertions 
which alter the bottomonium energy levels. The loops are composed of a charm and anti-charm 
pair, the shaded ellipses represent cumulative second order zero mass quark and gluon effects, and 
the double lines depict either a b or b quark. 



functions in one dimensional integral representation will be the basis for our calculations 
later. In section 3, we shall briefly display the Green's functions for arbitrary angular 
momentum, but the complete treatment is devoted to the S state, as we shall use it in a 
subsequent section. Our momentum wave functions will be used in section 4A to obtain 
the ground state energy level shifts for IS bottomonium due to next leading and next to 



next leading order charm mass effects in single insertions. We are able to give our results 
completely analytically in simple functions. They agree with what were obtained partially 
analytically by others. In section 4B, we give our formulation for double insertions and 
carefully state our subtraction procedure. A simple test is to apply it to obtain the ground 
state energy level of the modified Coulomb potential — ~a s (l + 5)-, the exact result of which 
is known a priori. Numerical analysis is presented in section 5 and some concluding remarks 
are made in section 6. In an appendix, we further illustrate our formalism by evaluating 
some well known matrix elements of (i) ' 1 ' 2 . 

II. MOMENTUM WAVE FUNCTIONS 

In this section, we shall derive a one-dimensional representation for the color Coulomb 
wave functions and show that, when converted back into coordinate space, they are the same 
as in textbooks. In fact, they will be slightly more general, because in some cases where the 
perturbing potential has power law dependence on the radial coordinate, we can take care 
of it easily. We shall use an algebra to determine the normalization factors 

From the Schrodinger equation, with the reduced mass m = rrib/2 for bottomonium, 

2m i <9£ o Pi 

we take away the centrifugal barrier and perform an angular momentum decomposition 

m = yuo.<pmoc {1+1 \ (2) 

to obtain 

'( £ -i(7^-^' + 1 O + 5 aAKH0 - (3) 

We have used £ to denote coordinates, with £ as the radial distance and 6, <fi as the angles. 
Now we perform a Fourier transform with respect to the radial distance £ 



to arrive at 



Q 



^)= dpe^jiip) (4) 

V 2tc J-oo 



Let us denote the energy of the state as 

2m' 
and rearrange the equation into 



E = A (6) 



where we have written 



d + i+i + s, + i+i^, = (7) 

dp p - p P + Po 



4 m 

^=o a ^ — • 8 

3 Po 

The solution of this equation is 

T'(P) = A (p-p )W+^(p + p )W-^' (9) 

where A^ is an integration constant to be determined shortly by normalization. The bound 
states are determined by the poles in the upper p-plane. Thus for the negative energy 
solutions we set 

Po — IK, i] = —in, (10) 

in which we require 

n -l-l = r = 0,1,2... (11) 

so that there is no pole in the lower p-plane. We do not want to burden our expressions in 
the remainder of this section with indices and therefore it is understood that we shall deal 
with states with the same principal quantum number n one at a time. 

To determine Ai, we go back to the coordinate representation, where we notice that 

«,(0 = r'7i(0 (12) 



satisfies the equation 



where 



Hm = Eu h (13) 



m p 2 l{l + l) A I Id 



Let us consider the operator jll] 



A =p-z(7- 3««— )• ( 15 ) 
Some straightforward algebra shows that we have 

LiHi - H^Li = 0. (16) 
It follows that if u\ is an eigenvector of H h so is L L ui of Hi_i, i.e. if the eigenenergy is E Q , 

Him = EqUi, (17) 

then 

UHiui = Hi-tLiUi, -> E (LiUi) = Hi^Lm). (18) 

£0 as well known depends only on n but not 1 for bound states. One can easily arrive at 

Therefore, for a given n, Li acts as a lowering operator 

Lm = B lUl _ 1} (20) 
which gives a norm 

5^1 ' = (f-)Vp--»' (21) 



/? 



upon using the bound state energies 



We shall follow the conventional choice to make the wave functions u\ real, which dictates 
the choice of the phase for Bi so that 



L1U1 = -i-\ /-j jUj_i, a x = -a.m. (23) 

a V I n z 3 

When we transcribe this last equation into the momentum representation, we have 



1/1 Id ,, 1, 

a^-^W l - 1 = itp -ra hl - (24) 
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Upon using the explicit solutions for 7(p)'s, we arrive at 



A t = vV^PAh, (25) 

na 



and upon iteration 



A t = (-2k)' 


(n + /)! 1 
n[n — / — !)! 


We have used the relation 

1 


K= — , 
na 


and in conformity with the standard choice of phase, we fix 


A = 





(26) 



(27) 



(28) 



At this point we want to be reminded that the radial wave functions are 7z(£)£~^ +1 )- 
Because we shall use them to evaluate matrix elements of operators which may have £ 
dependence, let us define more generally 



the Fourier transform of which 



~" l °-7?(fl, (29) 



7*(0 = _/ dpe«7f(p) (30) 



satisfies 



and therefore 



^7f(p)=7T 1 (p), (31) 



(^)*7f(p)=7?(p) = 'W(p)- (32) 



It can be inverted to yield 



-i ! /-° ju^t-iA (P' + P + Po) 11 - 1 - 1 



«'*) - (-d ,+1 (^t)j y> ^+ 17,™^ • (33) 

for a given n and £ > 1. This is the basic integral representation we shall use, although it 
can be evaluated into Gegenbauer polynomials. 



To proceed further, we want to write Laguerre polynomials L^ in an integral representa- 
tion over momentum. From the definition of 7(£) and how L k r appears in the wave function, 
i.e. 



tynlm J-lm\"i 



Y lm (0A) 



f 2 ^ 3l i^r e ^^^^ 2 ^' ™ 



where r = n — I — 1 and k — 21 + 1, we have 



(r + k)\ 1 g f° ,,. A _, 



4?(p - 2«e) = -^^ (^i)! e< / d p'p' 



X' 



2TTiJ-oo F {p + p' -J9 ) r+fc+1 ' V ' 

One can of course carry out the p and p' integrations to obtain 

L k M = ^P' k (^-/p r+ke ~ P - ( 36 ) 

However, for some problem where a sum over r has to be performed, it is more useful to 
have Vl presented as in Eq. (}35j) . 

Before we leave this section, we give a formula which is useful in the calculation of matrix 
elements, 

p p\p{ r(z + i)r(j + i) f°° p'^ +1 

/ dpi dp 2 t ■ —rrk = -rf | -| n / d P f , , , u , 37 

Jo (pi + p 2 + o) r(« + j + 2) Jo (p' + o)* 

when k > i + j + 2. In an appendix, we shall use it to reproduce some of the classic matrix 
elements to illustrate further our formalism. 



III. PROPAGATORS 

In order to carry out double insertions for the next-next leading order energy shifts, we 
need the propagators which are defined in 

G(x, y; E) = -Y: ^^ = £(2Z + l)G t (x, y; E)Pi{—), (38) 



where Pi are the Legendre polynomials and please note the minus sign in the sum over 
the eigen-energies E{. We have assumed a central (Coulomb) potential in writing down 
the above expression. Because of an interesting result in angular momentum decomposed 
Green's functions |l3j , we make a slight change of notation and write the energy as 

E = -«o/m 6 = Pl/mb, (39) 

then 



G,(xvE) - mfc/t0 f2:: x) l (2Knv) l e^^ V I>? +1 (2«ox)Lg t + 1 (2« y)r! 
G^x^E) - ^ (2k oX ) (2K y) e 2^ {r + l + l _ v){r + 2l + l y 

m b K ™ 2n 1 u l n (2Kpx) u l n {2n y) 
2n ^ 1 {2K fn-v x y ' [ } 

where we have restored the principal quantum number n to label the radial wave functions 
(i.e. u\ — > u l n ) and defined 

" s 2^ ( 3 a>) - (41) 

We shall later on be interested in the behavior v — > 1 when we look into the ground state. 
Now because of Eq. (|53|) , which gives the reduced radial functions 



m " (2k ° x) = (i)~(w)(-iyi^(-L) 1/2 

and a similar expression for u "^ Koy ' , with pi — > p 2 an d p^ — > p 2 - Also, we write 

1 r dpp"-"- 1 , (43) 



n — v Jo 
assuming 1 — v > 0. Performing the sum over n in Eq. (j4(J|) . we find 



oc 



G l (x,y;E) = (-l)'^^^^) 2 / "" dpie**« / (W™ 



4vr (/!) 2V 2vri 
x f dp' lP 'l f dp' 2 pl ^ dpp l - v 

J—oo J — oo JO 



oo J—oo 



X 



((Pi + Pi - Po)(P2 + P 2 - Po) - P(Pl + Pi + Po)(P2 + P 2 + Po))^ 1 ' 

fc = 21 + 1. (44) 
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For the case / = 0, it is easy to perform the sum and integrate over p'i,p' 2 to obtain 

rn, 1 roc roo 

G l=Q (x,y;E) = ^i,—) 2 / d Pl e^ x / dp 2 j™ 

xj o dpp-^lni-), (45) 



where 



and 



be = (pi - p )(p2 ~ Po) ~ p(2piP2 ~ 2p ) + p (pi + p )(P2 + Po), (46) 



ad = bc- AppQ. (47) 

With the wave functions and the propagators, we can proceed to calculate energy shifts 
in the next sections. 

IV. ENERGY SHIFTS DUE TO CHARM EFFECTS 

One objective in high precision tests of fundamental physics is to give operational meaning 
to the parameters in a theory and to measure them. Among the many parameters in the 
Standard Model, great improvements are being made in extracting the values of heavy 
particle masses, both experimentally and theoretically. For example, for the b-quark mass 
one way to determine it is by using the ground state of the bottomonium. As for the top 
quark, although its life time is too short for forming toponium, the location and the shape 
of the tt threshold will yield crucial information. 

Such considerations have been augmented to a very high and sophisticated degree by 
many groups. Thus, the ground state energy level of the bottomonium has been calculated 
to an accuracy of a^.(ln(a s )) [1]. At this order, among other contributions, there are the 
shifts due to the charm quark vacuum polarizations in the potential function. The important 
ratio here for the energy shifts is k/m c , where k is the momentum transfer, which is ~ |a s ra&. 
Thus k/m c ~ 1, in contradistinction to what is in atomic physics, where it is ~ a em . This 
numerical value is a cause for concern if one is to perform the intended calculation by 
an approximate expansion either in k/m c or m c /k. We must calculate effects due to these 
potential terms exactly. In this section, we shall show how this is handled in our formulation. 
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A. Potential Terras and Energy Shifts for Single Insertions 

The expression for the energy shift due to a single insertion is 

d 3 m&v(om, (48) 

-oo 

with the Fourier transform 

1 /•oo ^ ^ 

V ^) = V~2 d 3 ke* k <V(k 2 ). (49) 

Z7T J — oo 

After introducing the momentum wave function as in Eq. fjHOj) . 

V /7T J — oo 

into the spatial wave function 

because the potential is spherically symmetrical, we can immediately integrate over the solid 
angle Q to obtain 

AE=- dd d Px ^H(j>i)V{Z) dp 2 e^H(P2)- (52) 

Z7T JO ./-oo J-00 

At this point, we use the momentum representation of the potential in eq. (|4T)|) . We write 
k ■ £, = k^cosO and integrate over £ and the angles to obtain 

ae=-^ r d Pll \{ Pl ) r d P2l \{ P2 ) r kdk ^i+^+v ^ (53) 

111 1 J-00 J-00 Jo P l + P 2 ~ k 

Let us specialize and consider the energy shifts of the ground state, the wave function of 
which is (n=l) 



*°°M = ^^F < 54 > 



with 



/ 2k i 4 2 ,«N 

^0 = -\/ , «i = ^a s m reduced = -a s m b (55) 

V 7T 6 A 



for a bottomonium. We can easily use the residue theorem to carry out the integrations. 
We shall write as our standard form for any piece of the momentum space potential 

12 



V (k 2 ) = l 2 9(k 2 ), (56) 

taking into account of the factor of 2tt 2 in Eq.(|49|). then it is easy to show that the corre- 
sponding energy shift is 

- < - g{ e )>=N j <s m-—^^), (57) 

where the factor N is 

iV = Agl6«i. (58) 

The potential pieces to account for charm effects are due to: (a) the charm loop in the lowest 



order 



4 ~ ^ - 2 r r ^ /(*) L(-hn(—\ + ~" 

' k 2 + 4m 2 z 2 k? K 2 U{ m 2) 6' 

where 



4 9 r°° f(?\ 1 1 k 2 Fi 

Vffi° <*>) = (- r ,)r P a,-( / d z >U^ + F (-^«(^) + g)), (59) 



/(*) = Iv?=T(i + -L), (60 ) 

b) the iteration of the above and its combination with the zero mass quark and gluon effects 



vl L ° less (k 2 ) = (-5«-)£j^( - AM^) + «i), (6i) 

where 

A) = y^a - gTj?m, oi = — CU - -£r T Fni, (62) 

and /i is a subtraction scale, C^ = 3, Tp = 1/2 and rfy = 4, and (c) a s corrections to the 



oop. In this section, the mass symbol m refers to the mass of the charm quark. We 



charm 
have J5, 



yNNLO _ yNNLO(-i\ , yNNLO^) _L yNNLO(o\ (co\ 

charm charm \ ) ' charm \ I ' charm \ /' \ UJ / 

where 



C°(l) = ex [ Sl jT <fe/(*)^|^ + a 2 /°° dz 



k 2 + 4m 2 z 2 J\ k 2 + 4m 2 ,2 2 

1 A; 2 1 A; 2 1 
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with 

a 2 

- Cl = (-4a s /3)(^)(T F a s -), 

4:771 

01 = -2/3 , a 2 = -2 / 3oZn(— =-) + 2ai, a 3 = /3 , 

2 2 

777, 777. 

a 4 = -(5/3)/3 + fioln(— =-) - ai, a 5 = -(5/3)/? M— 7) + ( 5 / 3 ) a i> 



V charm \ l ) ~ C 2[l (iZ x dZ 2 A o„2~ L " 2 'J \ Z 



where 



1 k 2 + 4.77i 2 z\ k 2 + 4m 2 z\ J\ ' A; 2 + 4m 2 z 2 

+(-2fn2 + 5/3)jr ^^¥ + 4P H#V) + 5/3)2]l (65) 



4 2 

c 2 = (--a s )(7>a s — ) 2 , (66) 

O D7T 



and 



V2SP(3) = 4m»(1 + i) + 4M1 + J^) - (H^) - J£ - »«], (67) 



where 



a.,,o,76„, Zn^ , Zn% 



and 



1 £?1 TO 

A = exp(— - + — C 3 -/n2), c 2 = 0.470 ± 0.005, d 2 = 1.120 ± 0.010. (69) 

228 19 

We now apply Eq. (|57jl to obtain energy shifts. The more tedious task in the calculation 
is to integrate over the spectral density f(z). It turns out that all the integrals involved 
can be evaluated analytically into simple functions. Our method is first to assume that the 
parameter 

y = ^ (70) 

m 

is < 1, so that we can expand the integrand in powers of y/z. The integration over z can 
then be performed and the infinite series are then resummed. Since we have combinations of 
integrals which correspond to physical quantities and the analytic results for them, we can 
then unambiguously continue them into values of y > 1, which are relevant for investigating 
the toponium and the very light quark mass limit. 
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For the next to leading order, we have 



where 



Q 



(A£)™ m = {-\oc s )T F a s^-K X {y*I - (ln(2y) - H)), (71) 



1.2 11. 1 cos- l y,2 1 

/ = ( 1 ) -( 1 1) 

?/ 2 V Q } (1-y 2 ) 1 / 2 y V 2y 2 } 

+1^(4 +|). v^ 1 ( 72 ) 

2y 6 y z I 

which goes to 2/5 as y — ► 0, reproducing the well-known QED Uehling result |lfi| . By 
analytic continuation, we have 



r _±r A + Hi _ i fa(y + (y 2 -i) 1/2 ) ,2 j_ _ 

y 2 V 6 J (^-l) 1 / 2 y V 2y 2 ] 

+ w4 + l } - s£i (73) 

which will be useful when we consider tt threshold or high Z muonic atoms. The result here 
a^ees with what was obtained hefore Q. 

For the next- next leading order, we just list the analytic result for each piece. They are 



1 r°° ln(k 2 /Am 2 ) r°° r°° k 2 ln(k 2 /Am 2 ) 

2^ < Ji i{Z) k 2 + Am 2 z 2 > A dzt{z) h \k 2 + AK{) 2 k 2 + Am 2 z 2 

2 f°° j si f°° j, > k ' 2 ln ( k ') 



*«! dz!(:) L dk ' 



(2m) 3 A Jyj Jo (k 12 + y 2 ) 2 k' 2 + z 2 

N J^S^ dzf{z) Ty [ ll^ 2Hy) 

2l,: m*)+ l 1 



[z 2 — y 2 ) 2 z 2 — y 2 



N Wfl^W H+(,nW2y ^ + l) + l)G] ' 



(74) 



j /*00 I /*00 p 



oo J.2 

dk- 



2-k 2 A k 2 + 4m 2 z 2 " ' A ' Jo {k 2 + 4/tf) 2 k 2 + 4m 2 z 2 
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and 



"^(f) ft/Ml/"'"" 1 



(2mf y 2y J Ji Jy n (z 2 -y 2 ) 2i 2 z 2 - y 2 

N 7^Y 3 (-r^ J -y G ^ ^ 

(2m) 6 Ay ay 



<T L/ n 2 £)> = N Hdk- /n2( - 2 ' 



2vr 2 k 2 K m 2 ' Jo (P + 4/t 2 ) 2 



71 r 7T 2 1 



N Ij j[- + -ln\2y)-ln{2y)l (76) 



1 <T L/ n (4)> = n Hdk- /n( - 2 



2tt 2 k 2 m 2 ' Jo (k 2 + 4K 2 



N ihi M2y) ~ 1] ' (77) 



2tt 2 P l + 4m 2 d 2 j io (fc 2 + 4/t 2 ) 2 



4m 2 d 2 ' 



lQm 3 y 2 y d y + d ' 



1 1 « T 7T ,„„, 

< 77 >= ^7^' (79) 



2tt 2 A; 2 32k? 



I r 00 f 00 A; 2 

— <^ cfc^ dz 2 f( Zl )f(z 2 ) (jfc2 + ^ zlw + Am24) > 

/■oo />oo 1 2 2 z 2 

<l dz! J dz 2 f(z 1 )f(z 2 )(- 2 — ^){ - 2 , ]_ 2 _ 2 - L2 , jz^ )> 



i ji zf — z 2 k 2 + 4m 2 z 2 k 2 + 4m 2 z 2 ' 

n ,00 /•° ,„ /N „ /Nr/ n(i 1 N/ o 1 1 



N T7^-T dz! dz 2 f(z 1 )f(z 2 )[(y 2 — + -)(-y 



lQm 3 y J\ L Ji '' v v Z/Lvy ch/ 2 2 z\ — y 2 z\ — y 2 

d / 1 ^2 2 Z l 1 \i 

y i 9\ i 2 9 y 2 92 9 /J 

Oft/ 2 2i + 2 2 ^2 - I/ 2 ^i - I/ 2 ^ - £T 



The functions G, H, J, K which individually can be defined only for y < 1 are 
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oo 



G = I <lz f ^ Z 



i z 2 — y 2 



1 i n_ ? ,2U/2 i 

?I |1+ V )(1 - » " ,W| -6l (81) 



H = r iz ?mw*) 



l z 2 — y 2 



1\ 



2 [4(1 - y 2 ) 1/2 Ml + (1 " 1/ 2 ) 1/2 ) " 3Jn(2) + i 
8?r 2 

+4((1 - y 2 ) 1/2 Jn(l + (1 - y 2 ) 1 / 2 ) - J„(2))], (82) 



j= r^ # 



i z 2 — y 2 
n 1 r 3 



2 , 2L4 (l-^) 1/2 + ^(l-(l-^) 1/2 )]- (83) 



/■oo 1 

K = / dzidz 2 f{z 1 )f(z 2 ) - 



Z\ + z 2 z^- y 2 
1 19 9 13 1 

+ ,»-. fe )(i -, 2 )./ 2( __L__I___L ) 

^-^ 1/2 ^ + 237 + ^ < 84 > 



Using them, we have the following results for y < 1: 



1 2 , , , 5 11 , . , ,ln(2) 3/n(2) 1 



3- 



+52(1 - y yi\\ + - 3 )+ ^i - y r i/2 (^- 3 + -)+gi(i- y 2 ) 1/2 (-^- 5 - - 3 ),m 

V V ^y 6 y 2y° y 6 



-f ( J - ^ = 4 + A + <-h - A) + 9i(i - y 2 ) i/2 (-4 - 4) 

ay y 4 by z y° Ay 6 y a y a 

+ 9l (l-!/ 2 r 1/2 (~-i), (86) 

17 



(^h**) -£<*-*»> 



7 13 85 2 2 7 15 1 , 

V + V + 3%2 + ^ 5y^ + ^ ^ 1%8 ~ 1%6 + 4^2' 



where 



and 



+(1 - vr^M— - ^ + £ + £) + »?<£ + ^ - £), (87) 



7T 

S^i = sin~ {y) - - = -cos~ (y), (88) 



2 



£ 2 = /n(j/) sirC\y) - | Zri(l + (1 - y 2 ) 1 ' 2 ) 

= -ln{y)cos-\y) - ^ln{- + ((-) 2 - if' 2 ). (89) 

This completes our evaluation of the expectation values of various pieces of potentials into 
simple functions. Please note that g\,gi and (1 — y 2 ) 1 ^ 2 appear together in the correct 
combinations to allow us to analytically continue the energy shifts into real functions on the 
proper branch for y > 1. They are given by the substitutions 

cos-\y) ln(y + (y 2 - l) 1 ' 2 ) 

v < i -► — r^ — Tvm — v> 1 , ( 9 °) 



(1-^)1/2 *- (2/ 2 -l)V 

and 

Ml + ((1)2 _ i)V2) cos-H 1 ) 

— 1W iV - - « < 1 _> V_ y > 1 (91) 

((1)2-1)1/2 V- 1 (l-(i)f 2 V - L ^ 

We have a complete analytical result for the single insertion of the potential to account 
for the finite charm mass effects for all values of y Some of the integrals were not given in 
analytical forms by other authors and can be calculated numerically [6]. The values agree 
with ours within numerical accuracy, although we are somewhat unsure how the numerical 
integration programs decide on what branch to follow for y > 1. 

B. Potential Terms and Energy Shifts for Double Insertions 

In the next-next leading order, we have energy shifts due to double insertions of V^^ 
and V^^f iess of Eqs. (l5T?j) and (pTTj) . respectively. For a state with wave function ipE n=1 and 
perturbing potential V, the shift is given by 

&E doub ie insertion = - d 3 x d 3 y ip En=l (x)V(x)Lim E ^ En ^ 1 G(x, y; E)V(y)ip En=1 (y}- (92) 
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Note that the regularization applies to the Green's function only. We shall discuss this point 
later on when we compare results. We write 

1 

2~ z 



V(x) - ^ d%^ s V{kl), 



1 2n r°° 

lit 2 ix jo 
1 x 2tt r°° 



(-^-)— / dkMe lkx ~ e~ ikx )V{kl) 

2lX z IX JO 



/oo 
dkxhe^Vikl) (93) 

-co 



2lT 2 ' IX J-oo 

and a similar expression for V(y) with k\ — > ki- Likewise, we express 

II /-oo 

V^) f = YL(Q x )(-)(-) 1 / 2 / dpe^Kp), (94) 

X 2ix J-oo 

and for ip(y) with p — > p' . Upon using the addition theorem 

Pi(^-f) ~E^)^(^), (95) 

xy £i -f i m 

to carry out the angular integration of x and y, we arrive at 

fi.jr /"OO /'OO /"OO /"OO /"OO /"OO 

AE douUe insertion = ~ /0 9 -, 9 / cfe / */ / ^1^1 / dk 2 k 2 / dp / dp' 

(Z7T z 2j z JO JO J-oo J-oo J-oo J-oo 

Xe ifc 1 x e tfa tfe i^ e i^( p ^^y( A .2)y( jfe 2j Gj ( a . )y; £^ (gg) 

Now we specialize to the ground state given by Eqs.(|54|) and (|55|). We integrate over x 

and y, by using 

^(fci+P+Mx = , ( 97 ) 

Kl + P + Pi 

and Eg. (|45|) . which give 

AP - 2m fe «i 1 2 /■ !_„ 

<-^ I - J double insertion /„ a\o ^„ ./ / "fp 

(27T Z ) Z K 2m JO 

/OO /*00 /*00 /*00 /"OO /"OO 

d^i^i / dk 2 k 2 / dp / dp' / dpi / <ip2 
-oo J— oo J— oo J —oo J —oo J —oo 

x(^-) 2 (^~ ) 2 , -±- tt-^TT *"(^)y(*i a M*g) 
p — zfti p' — z«i ki + p + pi k 2 + p + p 2 ad 

2m b k\ f l 



rl /"OO /"oo /"oo /"oo 

/ dpp " / dk\k\ I dk 2 k 2 I dpi I dp 2 

JO J —oo J —oo J —oo J —oo 



(2n z ) z K JO J-oo 

1 1 he 

x (- — ) 2 (- — ) 2 M"^)^(*i)^(^)- (98) 

To organize our calculation better in what follows and so as not to have to track the 
branch cuts of products of logarithmic functions, we shall write 

TrAf r 0/72 . . 4 ._ 2 T . f°° , /(z) 1 r 1 dt 1 



1 1 m 2 5, 



+ Ii(oM-)+z)). (99) 
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and 



2 



Vl L ° less (k 2 ) = (-~a s )^Lim^ ( - ft [ J k 2 + ^ /t + ^ { ~^ H ^ ] + fll)) ' (1 ° 0) 
Therefore, it is clear that after we calculate the energy shift due to 

4 2 r°° 1 

V^k 2 ) = (--a s )(a s T F — )jf dzfiz)^^, z = 2mz } (101) 

we can make minor changes on the spectral density and take z 2 — ► fif/t or z — > to get the 
other contributions. This is what we are going to do immediately. 
The integration over k\ in Eq. (j98j) gives 

/oo I I I 
dk x k x -„ A-. — f = m{ -) 2 , (102) 

-oo kf + zf ki+pi + lKi Pl + l{Kl + Zi) 

and a similar result for k 2 . This leads to integrations over p\ and P2, which give a factor of 

(2m) 2 -^—ln(^Y (103) 

dp\dp 2 ad 

evaluated at p\ = —i{ni + Z\) and p 2 = — %{k\ + z 2 ). 

We are interested in the limit of 

p = Lim e _ +0 m 1 (l + e), e = 1 — v. (104) 

Thus in the calculation of the energy shift, we shall keep only the e° and the - terms. Since 
the p integration will produce the - pole, we must keep terms to order e 1 in the integrand 
as well. Now it is also clear that all terms of order ep n with n > 2 can be dropped. We 
shall show later on that the - terms cancel, because we are interested in the subtracted 
propagator 

Lim E ^ En=1 (-G{x,y;E) - \E n=1 > — < E n=1 \), (105) 

-& — &n=\ 

to obtain the regulated level shift. 

Carrying out the differentiation and evaluating at the value of pi )2 as discussed, we find 



d 2 , ,bc 



M-j) = p 



dpidp 2 ad 

< 



An 2 8k 2 

— e- 



((2«i + zi)(2«i + z 2 ) - p Zl z 2 ) 2 (2«i + z 1 ) 2 (2k 1 + z 2 ) 2 

8Kf(2/ci + ^i + 2Ki + z 2 ) 
(2«i + ^i) 3 (2 Kl + z 2 f 
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(106) 



We are now faced with an integral 



Ak\ f dpp- v 
Jo 



1 



+ e- 



((2«i + zi)(2«i + z 2 ) - P^ 2 ) 2 (2«i + ^i) 2 (2 Kl + z 2 ) 2 
2«i(2«i + Zi + 2«i + z 2 



— e- 



(107) 



(2«! + Zi) 3 (2«i + ^ 2 ) 3 

which can be easily calculated by adding and subtracting the terms in the square brackets 
at p = 0. Putting these pieces together, we obtain 

- 3 ',4 w „ 2 \ 2 

1 



A£(VTO = 2m 6 ^ (-a s )(a s T F — 



/OO /"OO 

dzifiz^ I dz 2 f(z 2 )4:K 2 1 



(2«! + %) 2 (2k 1 + ^ 2 ) 2 



_1 (2«i + zi)(2/ci + gg) ~ ^2 v 

e nl (2k 1 + z 1 )(2k 1 + z 2 ) j 



Z1Z2 



(2«i + ^i)(2ki + z 2 ) - ^i^2 
1 



+2«i 



+ 



(2 Kl + z 1 ) 2 (2fi: 1 + ^) 2 
1 \ 



(2«i + ^i) 2 (2k 1 + z 2 ) 3 (2«i + Zi) 3 (2«i + ^ 2 ) 2 / 
There remain two items to be settled. We first expand 



(108) 



K 



-1 .--1 



Ki\l-e) 



(109) 



in Eq. (jl08|) and then perform a subtraction as indicated in Eq. (jl05|) because the ground 
state should not be included in the propagator. This yields the subtraction term 

1 



d 3 x / d 3 y^{x) En ^Vi{x)^{x) En=1 



E — E n=1 



^{y)E n=1 Vi{y)^{y)E n ^ 



0(i)(l-| f )(A^- 2 



where 



4 9 r°° I 

{^E) NLO = (--a s )(T F a s -)(4K 3 ) / dzf(z)—- 



(2«i + *) 



^2" 



The end result is 



AE^Vt)^ 9 = 2m b KJ((^a s )(a s T F ^)) 2 

O 07T 



1 



X A ^ l/( ^A ^ /(Z2)4/C H (2«i + ^i) 2 (2k 1 + ^ 



x 



1 (2/tj + gi)(2Ki + 2:2) - Z 1 Z 2 , 

2 ™ l (2k 1 + z 1 )(2ki + z 2 ) ■ 



(110) 



111) 
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£i£2 

(2«i + zi)(2ki + z 2 ) - z\z 2 
1 



+2«i( 



+ 



(2«i + ^i) 2 (2ki + z 2 ) 2 
1 ,1 



(112) 



' (2«! + ^i) 2 (2/t 1 + z 2 y (2«i + 2i) 3 (2 Kl + z 2 ) 2/ J 
As a check that our regularization is correct, we look into the modified Coulomb potential 



4 ,1 .N 1 

•-a s (l + <>)-, 

6 X 



V(x) = --a s (l + 5)-, (113) 

3 x 

with the term proportional to 5 treated as a perturbation, which produces the ground energy 

E ground = -^a s ) 2 (l + 5) 2 . (114) 



4 v 3 ' v 

The <5 2 term is recovered from Eq. (|112j) by setting z 1)2 ~ ^ and requiring the spectral density 
to satisfy (a s T F -^) J dzf(z) -> 1. 

We shall take this as our standard subtraction, which will be applied to Eq. (|l()8jl and 
remaining energy shifts, by the operation 



H 



Kl, 



;ii5) 



To finish the calculation due to Eq.(|99|). we define 



V, 



t (A; 2 ) = Vl(fc 2 ) + K 2 (fc 2 ) + K 3 (fc 2 ), 



V 2 (k 2 ) = ( — a s )Tpa s — Lim u .^o( / — — T — % 

K J K 3 ' 3vr m> v 2 Jo tk 2 + fi 2 /t' 



4 
T 



2 
■3tt" 



1,1, ,m 



V 3 (k 2 ) = (-^a s )T F a s ^Lim^ -^Qln(- ? , , ,, 



+ 



and 



then the shift due to V\V 2 is 



bb = 2m b —(-a s (a s T F —)) 2 , 



l dt 



(116) 
(117) 

(118) 
(119) 



/■CO /■! WJ 

AE(2V 1 V 2 ) = 2(— 1/2)66 / <fei/(*i) / —4k? 

ji jo £ 



X 



+ Ln( 



(2K 1 +Z 1 ) 2 (2K 1 + fI i /Vi) 2 
(2/ti + Zi)(2/ti + /ij/A/t) - Z^j/y/t 
(2K l + Z l ){2K l +ll i /yft) 



+2 Kl 



(2«i + 2i)(2«i + /ii/v^) - zitn/Vi\ (2«i + ^i) 2 (2k 1 + ^/V*) 2 
1 1 V 



. (2«i + zi) 3 (2 Kl + A«i/>/t) 2 (2«i + ^i) 2 (2k! + /W*) 3 , 



(120) 
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We perform the t integration to obtain 



AE(2V 1 V 2 ) 



n d:f{z) \(2^Tw 



2Li 



_^_)_ 2 (i + Zn(-^)) 
2m + z J v 2 y 2 Kl JJ 



2 -{l + H^)) 



(2ki + z) 6 2k\ J 



121) 



Similarly, the V\V 3 shift is arrived at by taking the zero mass limit of z 2 , which leads to 



AE{2V 1 V 3 ) = bb r dzf(z)i—-^ 

Jl ZKi - 



-(-1-2) + —** 



,m 



(2 Kl + iF l( 2 fa( ^ ) + 6 ))l ^ <122) 



It is important to note that Zn(/ij)'s cancel in Eqs. (ll21lll22|) . as they must, because /ij is 
introduced to facilitate our calculation. Physical results should not depend on it. Shifts 
due to other pieces are deduced in a similar way upon being careful not to interchange the 
/ii — > limit and the t-integration. We have 



AE(V 2 V 2 + 2V 2 V 3 + V 3 V 3 



bb- 



44 



1/11 , . m . 



, m . 11 



7T 



- 6 - T + i 



(123) 



and 



/oo 



1 



m 



(2«! + zf 

4/ti 



2Li 



" 2in W - 3 ) + (2«x + z)3 



2 , 2 /ll m , 

•2«i + z' e V 6 v 2«i' 



, . m . 11\ 
2kx / 



(124) 



In a similar manner, we introduce 



vZssUn = vl(k') + v>(k 2 ), 



(125) 



where 



and 



and also 



Ku2\ 



V[{k 



V>(k 2 ) 



4 ,a s 

-a s ) — Lirrin-- 

3 Mvr Ml 


f , f di ! 


-«( ^y tv+rt/t 


4 a s 
--a s ) — Lim u . 

3 ; 4vr Ml 


- > o(^(-/?o/n(4) + a i): 



a& = 2m 6 ^(-o; s ) 2 (rFa s — )(-i). 

Kq 3 37T 47T 



(126) 



(127) 



128) 
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We have 



/oo 
dzf(z) 



20 o Li( 



2ki + z) 2 \ 2ki + z' 

~ (2/3q + a, + 2(3 Q ln{^)) - 2 ^H^) ~ A) - at 

(2ki + ^J- 3 V lK\ 



(129) 



and 



AE(2(V r 2 + \/ 3 )(K + ^)) = 2(o6)— j 



4k? 



-7(¥ + ' n O)( 2 * +a ' +2 ' 3 ° te( £; ) ) 



+ «(^) + ^ + t- % ( 6 -t + 1 , 



(130) 



We see again that ln(//j)'s cancel. The subtracted results of Eqs. ([12111130)1 as said are 
given by the prescription of Eq. ([115)1 . This finishes our double insertions. We have not 
been able to evaluate all the two dimensional integrals into simple functions as we could 
in section 4A. In fact, for y=l, some of these integrals produce the Catalan number, an 
indication that they are probably related to hyper-geometric functions. Our end result here 
agrees with that in ref.p. 

V. NUMERICAL RESULTS 

In this section, we present some numerical results. The mass in Eq.Q is the b quark 
pole mass M^ oe and the charm quark mass used throughout is MS mass m c (m c ). For our 
numerical work, we take 



m c = 1.25 GeV, M v h ° le = 5.0 GeV. (131) 

The typical energy scale involved is taken as \i = ^a s M^° e = 2.0 GeV, at which the running 
coupling constant is evaluated to be 



a a (2.0) = 0.30. 



(132) 



These give y = 0.80. The energy of IS state of bottomonium due to the non-zero charm 
quark mass is shifted by 
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FIG. 3: This figure compares the behaviors of AE /otg against AE /a? as functions of 

y = K\jm c . The solid line represents AE NNLO /af, while the dashed line represents AE NLO /a 2 s . 



A£f/° = -18.9 MeV, AE? s NLO 



-48 MeV, 



(133) 



due to next leading order and next to next leading order corrections, respectively. Altogether, 
these amount to a shift of bottom quark IS mass [17J. 



AMl s = -33 MeV. 



(134) 



With our explicit expressions for energy shift, we can calculate it for all values of y, 
including y > 1 . We display in Fig. 3 the behavior of AE NNLO /oP a and AE NLO /a 2 s as 
functions of y to succinctly summarize the intricate interplay of masses, by eliminating their 
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FIG. 4: The IS level shifts due to the charm mass effects are shown as functions of y at NLO 
(dashed line) and NNLO (dotted line). The total contribution of NLO and NNLO is represented 
by the solid line. 

dependence on the running coupling constant. In Fig. 4, we plot the NLO, NNLO and the 
total contribution to IS level energy shift as functions of y. 

Note that when m c — > 0, which corresponds to y —* oo, the energy shift due to double 
insertions tend to zero. This point was emphasized by Hoang p, because the leading order 
effects of the charm potential (Eg. ([59))) has no spatial dependence and the orthogonality of 
the wave functions which make up the propagator dictates this behavior. 

The ti "energy level shift" due to non-vanishing bottom quark mass can also be easily 



2(3 



estimated. We use 

M t = 174 GeV, m b = 4.2 GeV, a s {M t ) = 0.108, (135) 

which yield y = 2.98. Up to NNLO, the "IS energy level shift" due to the bottom quark 
mass effect reads 

AE = -25 MeV, (136) 

which is suppressed by the smallness of the coupling constant as compared with bottomo- 
nium. 

VI. CONCLUDING REMARKS 

One motivation for this work is to give a formulation to treat bound state and threshold 
physics completely in momentum space. For performing detailed high order calculations, it 
seems that the momentum space may have a natural setting, because Feynman rules are 
normally derived and given there. 

We have accomplished the construction of Coulomb wave functions for all quantum num- 
bers. Our propagators are quite compact, which among other things replace sums over 
principal quantum numbers by a parametric integral. The limit as the energy of the propa- 
gator^) approaches any of the bound state can be isolated easily. 

We have applied this formalism to investigate the charm effects on the IS bottomonium 
level shifts to the next-next leading order. Here, we have expressed all the single insertion 
results analytically in elementary functions. For double insertions, where the ground state 
is to be omitted as an intermediate state, it corresponds to a pole subtraction in -. 

Clearly, our formalism is useful for other problems. We have, for example, cursorily looked 
into the shifts due to charm effects in the next leading order on arbitrary bottomonium 
energy levels. We find the needed technique to be a slight extension of what we shall present 
in the appendix. Since the results for the shifts are known J8J , we shall not pursue it further. 
The S-states due to a linear potential can also be easily treated in this formulation. 
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APPENDIX A: MATRIX ELEMENTS IN QUANTUM MECHANICS 

We want to show how some matrix elements in quantum mechanics can be reproduced 
in our approach. Let us look at an integral for a fixed n 



/#f = f«H*) 



/ dp x pi' 1 / dp 2 p s 2 

J—oo J— oo 



A? 1 1 fO fO 

{i) s + t{ ' (s-i)\(t-i)\J- 

°° d (P1+P + Po) n ' l ' 1 (P2 ~P + Po)"'^ 1 , Al x 

-00 P (p! + p - p ) n+/+1 (P2-P- P ) n+l+1 ' 



Using the residue theorem, we have 



d4 e e ~ w s+t[ j ( s -i)!(t-i)!loo dpiPl 7-oo^ 2P2 

27TZ _j v ^ (2po + x)^- 1 (p 1 +p 2 -x)"-'- 1 _ 
(n + /)! l cfe J ( Pl + p 2 - 2p - x) n+l+1 lx -° ' l J 

If 2(/ + 1) > s + £, we can interchange the differentiation and the integrations over p 12 . 
Then, Eq. ()37jl gives 



fa 



tKOtKO = A * 2ni 1 ( d +l t-1 

e £* " (i) s+i (n + Z)!r(s + £) W 



/•oo 

x / dp' 
Jo 



(p> + &)"+'+i 

A| 2m "z^ 1 (-2p ) n ~'~ 1_i (" -1-1)! 
(z) s +*(n + /)! £- i!(n-Z- 1-i)! 

x y^ 1 ^ ( ^+i^+t-2j-H-2| a . _ q 

~^ j\(s + t — 1 — j)\n + / — i — j dx 

(A3) 

where b = 2po + x. Note that in most cases, only a small number of terms will survive the 
differentiation. For example, we can check the normalization Ai by setting s = t = I, which 
leads to 
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y l l > x CA4) 

p j\(2l-l-])\n + l-i-j (n-i + l)(n-i + l-l)...(n-i-l + iy v ; 



and then 



00 ,, 7* (0 7j(0 .. I" 00 ^„;,^„,/, 



/ #«(£) = / <%—# — tt= d Pimii{-p) 

JO JO C C J —oo 



M 2m «gi ( _ 2 ^ )n _ i _ 1 _.( n -/-l)!(n-i-0 



(i) 2 '(n + 0! S (n-i + Z)H! 

x(-^-) n+/ 6 i - 2 U =0 . (A5) 

ax 

Clearly only i = and i = 1 contribute to the sum because of the d/dx operation. We have 
after some simple algebra 



/■OO POO 

/ ^«(0 2 = / di 

JO JO 



7i (0 7i(0 



o v ^ ^o * e* e 



A ^(^f +3 7 -^ T (n-/-l)!(2n) 



(i) 2 ' v 2/€ y (n + 0! 

1, (A6) 



as it must. 



Next we consider s + 1 — 21 + 1, which will give us the expectation value of l/£. This 
time we use in Eq. flA3|) 



>p [ l > l - CA7) 

f^j\(2l-j)\n + l-i- j (n-i + l)(n-i + l-l)... (n-i-l)' 



which gives 



jf de«(o^(e) = / °°«^^ = £ ) *7f fl (p)7i(-P) 



A 2 2.i -- 1( _^ r _^( n - z -i), 



(i)«+i („ + /)! ^ ^ ™ (n-i + Z)!i! 

x(^-)™ + '6 l - 1 U =0 . (A8) 

ax 

There is only one term i=0 in the sum which is not annihilated by d/dx and we come up 

with 
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1 \2nJ n ~ l ~ i y-r l ^2^+2 



/ di u(Z)-u(Z) = A 2 2ir- 
Jo t 



= 4-. ( A9 ) 

n z a 
which agrees with a well-known result. 

As we stated earlier, to interchange the integration and the differentiation operations in 
Eq. (jA3|) . it is necessary that the condition 2(7 + 1) > s + t should be satisfied. For the case 
s + t = 21 + 2, which is relevant for the expectation value of l/£ 2 , we must treat it differently. 
Thus, we write 



f°, , m_m Af 2m 1 

k ? e +i e l+l ~ (i) 2i+2 (n + i)\(2i + i)\ ' K ' 



where 



We rewrite 



rOO Jt 7i(07i(|) = A 2 2ttc 1 

fi {iy+i{ n + l)\{2l 

B- /"^^(^^ ^- U (All) 



oo ,/ / ,,' „^n--I— 1 

o ^ ^ y dx' " {j)' + b) n+,+ ] 



p ' + x=p ' + b- 2p , (A12) 

and expand in binomial the numerator 

(p' + x) n - 1 - 1 1 (n-l- l)(-2p ) 

(p' + 6)™+'+! ~ (p> + 6) 2 '+ 2 + (j/ + 6) 2 '+ 3 

(n-/-l)(n-/-2)(-2» ) 2 , A , 

+ ~ WTVj^ + -- (A13) 

Only the first term will be non- vanishing after the d/dx differentiations, since all the others 
terms are convergent enough that we can interchange the integration and differentiation 
operations and are proportional to (d/dx) n+l b k , k < n + I. Then effectively 

After examining a few values of n and 1, one can easily convince oneself that 

B = -(n-l-l)\(2l)\b- {2l+1 \ (A15) 

and thus 
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o * ^ +1 e +1 nW + 1/2' { ' 



which again is a well-known result. 
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